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Abstract 

The production rates and substructure of jets have been studied in charged cur- 
rent deep inelastic e + p scattering for Q 2 > 200 GeV 2 with the ZEUS detector at 
HERA using an integrated luminosity of 110.5 pb _1 . Inclusive jet cross sections 
are presented for jets with transverse energies > 14 GeV and pseudorapidi- 
ties in the range — 1 < rf ct < 2. Dijet cross sections are presented for events 
with a jet having > 14 GeV and a second jet having E 1 ^ > 5 GeV. Mea- 
surements of the mean subjet multiplicity, (n s bj), of the inclusive jet sample are 
presented. Predictions based on parton-shower Monte Carlo models and next-to- 
leading-order QCD calculations are compared to the measurements. The value 
of a s (Mz), determined from (n s bj) at y cnt = 10~ 2 for jets with 25 < < 119 
GeV, is a s (M z ) = 0.1202 ± 0.0052 (stat.) t° mi9 (syst.) ^oo! ( th -)- Th e 
mean subjet multiplicity as a function of Q 2 is found to be consistent with that 
measured in NC DIS. 
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1 Introduction 



Measurements of the charged current (CC) deep inelastic scattering (DIS) cross section 
at HERA [UiaElllj at high virtuality, Q 2 , of the exchanged boson have demonstrated 
the presence of a space-like propagator with a finite mass, consistent with that of the 
W boson. Jet production in CC DIS provides a testing ground for QCD as well as the 
electroweak sector of the Standard Model. Up to leading order in the strong coupling 
constant, a s , jet production in CC DIS proceeds via the QCD-Compton (Wq — > q'g) 
and W-gluon-fusion (Wg — > qq') processes in addition to the pure electroweak process 
(Wq -> q'). 

At HERA, multijet structure has been observed in CC DIS 0E] at large Q 2 and jet 
substructure has been studied using the differential and integrated jet shapes 0. Another 
useful representation of the internal jet structure is the subjet multiplicity [HE]. The 
lowest-order non-trivial contribution to the subjet multiplicity is of order a s , so that 
measurements of the subjet multiplicity provide a direct test of QCD. 

This paper reports a detailed study of the hadronic final state in CC e + p DIS. Differ- 
ential cross sections are presented for both inclusive jet and dijet production. The jets 
were identified in the laboratory frame using the longitudinally invariant kx cluster al- 
gorithm After describing experimental conditions and the theoretical calculations, in 
Section 17.21 the inclusive jet cross sections are presented as a function of the virtuality 
of the exchanged boson, the jet pseudorapidity, r] :]ct , and the jet transverse energy, E^. 
In Section 17. 3| the dependence of the dijet cross sections on Q 2 and the invariant mass, 
m 12 , of the two highest-E^ jets are given. In Section \7Al the mean subjet multiplicity, 
(n S bj), as a function of the resolution scale, y cut and E^ using the inclusive jet sample 
is presented. Parton-shower Monte Carlo (MC) calculations and next-to-leading-order 
(NLO) QCD predictions fU] are compared to the measurements. In Section |H1 the value 
of a s (Mz) determined using the measurements of (n s bj) as a function of E^ is given. In 
Sectional the measurements of (n s bj) as a function of Ej, and Q 2 are compared to the 
results obtained by ZEUS in neutral current (NC) DIS [TT] . 



2 Experimental conditions 

The data sample used in this analysis was collected with the ZEUS detector at HERA 
and corresponds to an integrated luminosity of 110.5 pb _1 . During the 1995-1997 (1999- 
2000) running period, HERA operated with protons of energy E p = 820 GeV (920 GeV) 
and positrons of energy E e = 27.5 GeV, yielding a centre-of-mass energy of 300 GeV 
(318 GeV). A detailed description of the ZEUS detector can be found elsewhere jT2]. A 
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brief outline of the components that are most relevant for this analysis is given below. 
Charged particles are tracked in the central tracking detector (CTD) which operates 
in a magnetic field of 1.43 T provided by a thin superconducting solenoid. The CTD 
consists of 72 cylindrical drift chamber layers, organized in nine superlayers covering the 
polar-angle 1 region 15° < 6 < 164°. The transverse-momentum resolution for full-length 
tracks is a(p T )/pT = 0.0058p r © 0.0065 ©0.0014/p T , with p T in GeV. The high-resolution 
uranium-scintillator calorimeter (CAL) jHj consists of three parts: the forward (FCAL), 
the barrel (BCAL) and the rear (RCAL) calorimeters. Each part is subdivided trans- 
versely into towers and longitudinally into one electromagnetic section (EMC) and either 
one (in RCAL) or two (in BCAL and FCAL) hadronic sections (HAC). The smallest sub- 
division of the calorimeter is called a cell. The CAL energy resolutions, as measured under 
test-beam conditions, are a(E)/E = 0.18 /\f~E for electrons and a(E)/E = 0.35/v^ for 
hadrons, with E in GeV. Jet energies were corrected for the energy lost in inactive mate- 
rial, typically about one radiation length, in front of the CAL. The effects of the uranium 
noise were minimised by discarding cells in the electromagnetic or hadronic sections if 
they had energy deposits of less than 60 MeV or 110 MeV, respectively. A three-level 
trigger OES] was used to select events online. 

The luminosity was measured using the Bethe-Heitler reaction e + p — > e + p , ~f. The resulting 
small-angle energetic photons were measured by the luminosity monitor j!6j . a lead- 
scintillator calorimeter placed in the HERA tunnel at Z — —107 m. 



3 Data selection and jet search 

The selection of charged current events for the present study is very similar to those 
described in detail in previous ZEUS publications [3] . The efficiency of the selection cuts 
is typically above 90% and the remaining backgrounds are negligible. 

The principal signature of a CC DIS event at HERA is the presence of a large missing 
transverse momentum, p'xi arising from the energetic final-state neutrino which escapes 
detection. The quantity p'r was calculated from 

i>T 2 = Px + Py = I ^ E% sin &i cos <t>i I + I Ei sin 0* sin fa 



The ZEUS coordinate system is a right-handed Cartesian system, with the Z axis pointing in the 
proton beam direction, referred to as the "forward direction" , and the X axis pointing left towards 
the centre of HERA. The coordinate origin is at the nominal interaction point. The pseudorapidity 
is defined as rj = — In (tan |), where the polar angle, 6, is measured with respect to the proton beam 
direction. 
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where the sums run over all CAL cells, i, E^ is the energy deposit and 6*j, fa are the 
polar and azimuthal angles of the cell as viewed from the interaction vertex. The total 
transverse energy, E?, is given by E? = ^Ei sin Q^. 

The inelasticity, y, was reconstructed using the Jacquet-Blondel method [TJj and corrected 
for detector effects as described previously 0. The detector simulation was used to derive 
corrected values ^T.cor and y COT . The corrected value of Q 2 , Q 2 ov1 was calculated in terms 
of j4 TiCor and y COT using the relation Q 2 COT = ^ cor /(l - y COT ). 

The following requirements were imposed on the data sample: 

• fa > 11 GeV and Q 2 col . > 200 GeV 2 , to ensure high trigger efficiency; 

• Vcot < 0.9, to avoid the degradation of the resolution in Q 2 near y ~ 1; 

• $t/Et > 0.5, to reject photoproduction and beam-gas background. For the dijet 
sample, this cut was reduced to p'r/Ex > 0.3 with the further requirement that the 
difference between the azimuthal angle of the missing transverse momentum and that 
of the closest jet was greater than 1 rad. This cut removed poorly reconstructed 
back-to-back dijet photoproduction events; 

• a vertex position reconstructed with the CTD in the range — 50 < Z < 50 cm, consis- 
tent with an ep interaction; 

• the difference, A0, between the azimuthal angle of the net transverse momentum as 
measured by the tracks associated with the vertex and that measured from the CAL 
be less than 1 rad. This requirement removed random coincidences of cosmic rays 
with ep interactions; 

• p^ ack /^ T > 0.1, where p^ ack is the net transverse momentum of the tracks associ- 
ated with the vertex. This condition was not applied if p'r > 25 GeV. This cut 
rejected events with additional energy deposits in the CAL not related to ep inter- 
actions (mainly cosmic rays) and beam-related background in which fa has a small 
polar angle; 

• the event was removed from the sample if there was an isolated positron candidate 
with energy above 10 GeV, to reject NC DIS events; 

• a pattern-recognition algorithm based on the topology of the calorimeter energy dis- 
tribution and the signals detected in the muon chambers was applied to reject cosmic 
rays and beam-halo muons. 

The longitudinally invariant kx cluster algorithm jH] was used in the inclusive mode fH] 
to reconstruct jets in the hadronic final state both in data and in MC simulated events 
(see Section EJ). In data and MC, the algorithm was applied to the energy deposits in the 
CAL cells and in the MC it was also applied to the final-state hadrons. The jet search 
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was performed in the r\ — <fi plane of the laboratory, starting with the CAL cells or hadrons 
as initial objects. In the following discussion, Eta denotes the transverse energy, 77, the 
pseudorapidity and fa the azimuthal angle of object % in the laboratory frame. For each 
pair of objects, the quantity 

dij = [(r]i ~ Vj) 2 + (fa ~ 4>j) 2 ] ■ min(E Tji , E T ,j) 2 

was calculated. For each individual object, the quantity dj = {Eta) 2 was also calculated. 
If, of all the values {dy,<ij}, d^i was the smallest, then objects k and / were combined 
into a single new object. If, however, dk was the smallest, then object k was considered 
a jet and was excluded from further clustering. The procedure was repeated until all 
objects were assigned to jets. The jet variables were defined according to the Snowmass 
convention [TU] : 

i i T i T 

where the sums run over all objects associated with the given jet. This prescription was 
also used to determine these variables for the subjets. For jets constructed from CAL 
cells, jet energies were corrected for all energy-loss effects, principally in inactive material 
of typically about one radiation length, in front of the CAL (see Section 0J). 

For the inclusive jet sample, all jets with E^ > 14 GeV and — 1 < if et < 2 were retained. 
For the dijet sample, at least one additional jet with Ej, > 5 GeV and — 1 < r/ et < 2 
was required. The upper rapidity requirement is made so that the jet is within the CTD 
acceptance for efficient background rejection. There are very few events with jets with 
sufficient E^ below the lower rapidity requirement. 

With the above criteria, 1865 events with at least one jet and 282 dijet events were 
identified. 



3.1 Definition of subjet multiplicity 

Subjets were resolved within a jet by considering all objects associated with the jet and by 
repeating the application of the A>r cluster algorithm described above, until for every pair 

of objects % and j, the quantity d^ was greater than d cut = y C ut ' (J$t) ■ All remaining 
objects were called subjets. The jet structure depends upon the value chosen for the 
resolution parameter y cut . For each sample studied, the mean subjet multiplicity, (n sh j), 
is defined as the average number of subjets contained in a jet at a given value of y cnt '- 

(n sh} (y C ut)} = ^— KbiiVcut), 
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where n* b j is the number of subjets in jet i and iVj ets is the total number of jets in the 
sample. The mean subjet multiplicity of the inclusive jet sample was measured for y cut 
values in the range 5 • 10 -4 to 0.1. The y CVLt range was chosen to be small enough to have 
mean subjet multiplicities larger than unity and large enough to avoid the degradation in 
resolution caused by the finite size of the CAL cells. 



4 Monte Carlo simulation 



Samples of events were generated to determine the response of the detector to jets of 
hadrons and to evaluate the correction factors necessary to obtain the hadron-level jet 
cross sections and subjet multiplicities. The CC DIS events were generated using the 
LEPTO 6.5 program [20] interfaced to HERACLES 4.6.1 [21] via DJANGOH 1.1 |22|. 
The HERACLES program includes first-order electroweak radiative corrections. The 
CTEQ4D [23] NLO proton parton distribution functions (PDF) were used. The QCD 
radiation was modelled with the colour-dipole model [21] by using the ARIADNE 4.08 
program [23] including the boson-gluon-fusion process. As an alternative, samples of 
events were generated using the LEPTO model which is based on first-order QCD matrix 
elements and parton showers. For the generation of the LEPTO samples, the option for 
soft-colour interactions was switched off since its inclusion results in an increase both in 
particle multiplicity and energy per unit of rapidity that disagrees with the measurements 
in NC DIS at HERA [2B] . In both cases, fragmentation into hadrons was performed using 
the Lund string model [2Z] as implemented in JETSET 7.4 [2H]- To calculate the accep- 
tances and to estimate hadronisation effects, the generated events were passed through 
the GEANT 3.13-based j2Hl simulation of the ZEUS detector and trigger. They were 
reconstructed and analysed by the same program chain as used for data. For both the 
ARIADNE and LEPTO event samples, a good description of the measured distributions 
for the kinematic and jet variables was obtained [30J. 

To correct the data to hadron level, multiplicative correction factors, defined as the ratio of 
the measured quantities for jets of hadrons over the same quantity for jets at detector level, 
were estimated by using the ARIADNE and LEPTO models. Parton-level predictions 
were also obtained by applying the jet algorithm to the MC-generated partons. These 
predictions were used to correct the NLO QCD calculations to hadron level (Section |3J). 

HERACLES 4.6.2 j2I] was used to correct the measured cross sections to the electroweak 
Born level evaluated using the electromagnetic coupling constant a = 1/137.03599, the 
Fermi coupling constant Gf = 1.16639 ■ 10~ 5 GeV~ 2 and the mass of the Z boson Mz = 
91.1882 GeV [3^] to determine the electroweak parameters. 
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5 NLO QCD calculations 



The NLO QCD calculations were obtained from the program MEPJET which em- 
ploys the phase-space slicing method (221 • This is the only available program providing 
NLO calculations for jet production in charged current deep inelastic scattering. The 
calculations were performed in the MS renormalisation and factorisation schemes. The 
number of flavours was set to five and the renormalisation (/j,r) and factorisation (hf) 
scales were chosen to be /1r = /ip = Q. The calculations were performed using the 
CTEQ4M [23] parametrisations of the proton PDFs, which are based on the MS scheme. 
The jet algorithm described in Section [3] was also applied to the partons in the events 
generated by MEPJET in order to compute the jet cross section and the predictions for 
the subjet multiplicities. The cross sections were evaluated using the same values for a, 
Gf and Mz as in the electroweak Born level of the measured cross sections (Section HJ. 
In addition, the mass of the W boson was fixed to 80.4603 GeV. 

Since the measurements correspond to jets of hadrons whereas the NLO QCD calculations 
correspond to jet of partons, the predictions were corrected to the hadron level using the 
MC simulations. The multiplicative correction factor (Chad) is defined as the ratio of 
either the cross sections or the mean subjet multiplicities for jets of hadrons to the same 
quantity for jets of partons, estimated using the MC programs described in Section 0] 
The ratios obtained with the ARIADNE and LEPTO models were in good agreement 
and the mean was taken as the value of Chad- The value of Chad is ~ 1-03 (~ 1.10) for 
the inclusive jet (dijet) cross sections. For the mean subjet multiplicity, Chad is 2.13 at 
Vent = 5 ■ 10~ 4 and 14 < < 17 GeV and approaches unity as y cnt and increase. 

The theoretical predictions were redetermined after changing the parameters as described 
below. In each case the difference between the redetermination and the nominal prediction 
was taken to be the uncertainty in the calculation associated with the parameter under 
consideration. 

• Proton PDFs: the CTEQ5M [33] and MRST [33] sets, rather than CTEQ4M [23] , 
were used. Also, a set of the MRST PDFs with a larger d/u quark ratio at large 
Bjorken x was used. The uncertainty in the cross sections was less than ~ 4% for the 
inclusive jet cross section, except for high E^ , where it reaches ~ 20%. It was less 
than ~ 10% for the dijet cross sections. The uncertainty was negligible for the subjet 
multiplicities; 

• a s {Mz): the a s (M z ) values of 0.113 and 0.119, corresponding to the proton PDFs 
CTEQ4A2 and CTEQ4A4, were used. The uncertainty in the cross sections was 
typically ~ 2%; for the mean subjet multiplicity the uncertainty was ~ 1% ; 

• hr: in order to estimate the effects of the terms beyond NLO, the scale fiR was varied 
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between Q/2 and 2Q, while keeping \ip fixed at Q. The uncertainty of the cross 
sections was less than 5%. The uncertainty on the mean subjet multiplicity was ~ 3% 
for y cut = 10~ 2 ; 

• Chad^ the hadronisation correction, Chad, was varied by half of the difference between 
those evaluated using ARIADNE and LEPTO. The uncertainty typically amounted to 
less than 1%(3%) for inclusive jet (dijet) cross sections. For the subjet multiplicities, 
the uncertainty was less than 3% for y cut = 1CT 2 . 

• s m i n : the cut-off parameter s min in the phase-space slicing was changed from the 
default value of 0.1 GeV 2 to 0.01 GeV 2 . This uncertainty was less than 1% in all the 
calculations and was neglected in the estimation of the total theoretical uncertainty. 

The total theoretical uncertainty was obtained by adding in quadrature the individual 
uncertainties listed above and is shown as the hatched band in the figures. 

6 Experimental systematic uncertainties 

A study of the sources contributing to the systematic uncertainties of the measurements 
was carried out. The following sources were considered: 

• the uncertainty on the absolute energy scale of the jets was taken to be ±1% for 

> 10 GeV and ±3% for lower Ej? values [33]. The resulting uncertainty was less 
than 5% (12 %) for the inclusive jet (dijet) cross sections and less than 2% for the 
mean subjet multiplicity; 

• the uncertainty in the reconstruction of the kinematic variables due to that in the ab- 
solute energy scale of the CAL was estimated by varying the energy variables measured 
with the CAL by ±3%. The uncertainty was less than 5% for all distributions; 

• the differences in the results obtained by using ARIADNE or LEPTO to correct the 
data for detector effects were taken as systematic uncertainties; they were typically 
smaller than 5% for the cross sections and smaller than 2% for the mean subjet mul- 
tiplicities; 

• the selection cut of $t > 11 GeV was changed to 10 GeV and 12 GeV. This gave 
a variation of the cross sections (subjet multiplicities) of less than 5% (2%). The 
uncertainty evaluated from the variation of other selection cuts was typically less than 
2%. 

For the jet cross sections, the systematic uncertainties not associated with the absolute 
energy scale of the jets and the CAL are not point-to-point correlated and were added 
in quadrature to the statistical errors. They are shown as the bars in the figures. The 
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uncertainty due to the absolute energy scale is point-to-point correlated and is shown 
separately as a shaded band in each figure. For the subjet multiplicities all the systematic 
uncertainties are point-to-point correlated and were added in quadrature to the statistical 
errors. They are shown as the bars in the figures. 

In addition, there is an overall normalisation uncertainty of 2.0% from the luminosity 
determination, which is not included in the results presented in the figures and the tables 
of the cross sections. 

7 Results 

7.1 Data-combination method 

Due to the different centre-of-mass energy of the two data sets used in the analysis, the 
measured jet cross sections based on each set are presented separately in Tables [T] to 

The measured jet cross sections, a^, were combined using the following formula: 

comb _ ^"300 • £300 + Q"318 " ^318 

0-318 ~~ 7 <rr th /rr th } 4- C ' 

-L-300 • 1^300 /°318 ) + ^318 

where is the luminosity and cr^ is the predicted cross section. The ratio 0300/0318 was 
obtained using the program MEPJET. The ratio obtained by ARIADNE is within 1% of 
that obtained by MEPJET and was the same, within statistical errors, as that measured 
in the data. All the systematic errors have been assumed to be correlated between the 
measurements. 

No dependence of the subjet multiplicities with the centre-of-mass energy was seen either 
in the data or in the theoretical predictions; thus the subjet multiplicities were calculated 
using the inclusive jet sample of both data sets. The measured subjet multiplicities are 
presented in Tables E] to 

7.2 Inclusive jet differential cross sections 

The differential inclusive jet cross sections were measured in the kinematic region defined 
by Q 2 > 200 GeV 2 and y < 0.9. These cross sections include every jet of hadrons in the 
event with E 3 ^ 1 > 14 GeV and —1 < 7^ et < 2. The differential inclusive jet cross sections 
as a function of Q 2 , r/ jct and are shown in Figs. □ to H Both the ARIADNE MC 
model and the NLO QCD calculation MEPJET give a good description of the measured 
inclusive jet cross sections. 
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7.3 Dijet differential cross sections 

The differential dijet cross sections were measured in the kinematic region defined by 
Q 2 > 200 GeV 2 and y < 0.9. These cross sections refer to the two jets of hadrons 
with highest transverse energy in the event with E^' 1 > 14 GeV, E^' 2 > 5 GeV and 
— 1 < if 6t < 2. The differential dijet cross sections as a function of Q 2 and the invariant 
mass of the two highest-.&r jets, m^, are presented in Figs. |U and E3 The NLO QCD 
calculation gives a good description of the measured dijet cross sections. 

7.4 Subjet multiplicities 

The mean subjet multiplicity, (n s bj), was determined using the inclusive sample of jets 
in the kinematic region defined by Q 2 > 200 GeV 2 and y < 0.9. The (n sbj ) values 
were obtained using every jet of hadrons in the event with E^ > 14 GeV and —1 < 
77 jet < 2. The results are shown in Fig. El as a function of y cut for different E? regions. 
In the region of small y cnt , the ARIADNE MC model gives a better description of the 
multiplicity than the NLO QCD calculation MEPJET. At larger values of y cut both the 
MC model and the NLO QCD calculation give a good description of the measurement. 
In the region of small y cnt values, fixed-order QCD calculations are affected by large 
uncertainties and a resummation of terms enhanced by ln(y cut ) jS] would be required for 
a precise comparison with the data. At relatively large values of y C ut, a NLO fixed-order 
calculation is expected [H] to be a good approximation to such a resummed calculation. 

The measured (n s bj) at y cut = 10~ 2 as a function of E^ is shown in Fig. [7] The measured 
mean subjet multiplicity decreases as Ej, increases. The overall description of the data 
by the NLO QCD calculations is good. 

8 Measurement of a s 

The sensitivity of the subjet multiplicity to the value of a s (Mz) is illustrated in Fig. 
which compares the measured (n s bj) at y cnt = 10~ 2 function of E^ with NLO QCD 
calculations obtained with different values of a s (Mz)- Both the measurements and the 
NLO QCD predictions of the subjet multiplicities have smaller uncertainties compared to 
those of the jet cross sections. Therefore, the measured (^sbj), rather than the jet cross 
sections, was used to determine a s (M z ) using the following procedure: 

• the NLO QCD calculations of (n s bj) were performed for the five sets of the CTEQ4 
"A-series" PDFs, which differ in the assumed value of a s (Mz)- The value of a s [Mz) 
used in each calculation was that associated with the corresponding set of PDFs; 
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• for each bin i in , the NLO QCD calculations, corrected for hadronisation effects, 
were used to parametrise the a s (M z ) dependence of (n sh j) according to 

[(n sbj ) (a s (M z )) ].=l+d a s (M z ) + C\ a s (M z ) 2 . (1) 

The coefficients C\ and C\ were determined by performing a x 2 -fit to the NLO QCD 
predictions. This simple parametrisation gives a good description of the a s (M z ) de- 
pendence of (n s bj) over the entire range spanned by the CTEQ4 "A-series"; 

• this parametrisation was used to extract a value of a s (M z ) in each bin; 

• in addition, a combined value of a s (M z ) was determined by a x 2_ fit °f Eq. (JIJ to the 
measured (n s ^j) values for all bins. 

This procedure correctly handles the complete a s dependence of the calculations (the 
explicit dependence coming from the partonic cross sections as well as the implicit depen- 
dence coming from the PDFs) in the fit, while preserving the correlation between a s and 
the PDFs. 

The uncertainty in the extracted values of a s (M z ) due to the experimental systematic 
uncertainties was evaluated by repeating the above analysis for each systematic check. 
The largest contribution to the experimental uncertainty was that due to the simulation 
of the hadronic final state. 

The theoretical uncertainties arising from terms beyond NLO and uncertainties in the 
hadronisation correction were evaluated as described in Section 03 These resulted in 
uncertainties in a s (M z ) of Aa s (M z ) =t° ^lt and Aa s (M z ) = ±0.0014, respectively. 
The total theoretical uncertainty was obtained by adding these in quadrature. Other 
uncertainties described in Section El were small and were neglected. As a cross check, a 
linear parametrisation of the a s (M z ) dependence of (n S bj) was considered; the change in 
the extracted value of a s (M z ) was negligible. 

The values of a s (M z ) obtained from the measurement of (n s ^) for various Ej, regions 
are in good agreement. The value of a s (M z ) obtained from the measurements of (n s bj) at 
Ucut = 10~ 2 for 25 < Ej^ < 119 GeV, a region in which the parton-to-hadron correction 
was less than 10%, is 

a s (M z ) = 0.1202 ± 0.0052 (stat.) IS (syst.) t° S (th.). 

This result is consistent with other recent determinations using measurements in NC DIS 
of inclusive jet (HZIISH] and exclusive dijet cross sections [HHl as well as measurements of 
(n sb j) HU and with the PDG value, a s (M z ) = 0.1172 ± 0.0020 gO]. 
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9 Comparison of subjet multiplicities in CC and NC 



The present measurements of subjet multiplicities in CC interactions are compared with 
the corresponding measurements in NC DIS The NC data were reanalysed in the 
same kinematic region as that of the CC analysis. 

The measurements of (n S bj) at the value of y cut = 10 -2 as a function of E^ in CC and 
NC DIS are compared in Fig. |H^i. The value of (n s bj) is slightly larger for jets in NC DIS 
than for CC DIS for a given jet transverse energy. The NLO QCD predictions behave in 
the same way as the data. 

The subprocess population and the phase space available for QCD radiation depend on 
Q 2 . The measurements of (n s bj) at y cut = 10~ 2 as a function of Q 2 in CC and NC DIS 
are compared in Fig. IHJd. The values of (n s bj) in CC and NC DIS are similar and are in 
agreement with the NLO predictions. The differences observed in the subjet multiplicity 
as a function of E^ can be attributed to the different Q 2 distributions of the CC and NC 
processes. 



10 Summary 

Measurements of differential cross sections for inclusive jet and dijet production in charged 
current deep inelastic e + p scattering have been performed and are corrected to the elec- 
troweak Born level. The internal structure of the inclusive jet sample has been studied 
in terms of the mean subjet multiplicity. The results are given for jets of hadrons iden- 
tified with the longitudinally invariant kx cluster algorithm in the laboratory frame in 
the kinematic region defined by Q 2 > 200 GeV 2 and y < 0.9. Inclusive jet cross sections 
are presented for jets with transverse energies E 3 £ > 14 GeV and pseudorapidities in 
the range —1 < rf et < 2. Dijet cross sections are presented for events with a jet having 
E 3 ^ > 14 GeV and a second jet having E^ > 5 GeV. 

The predictions of the ARIADNE MC model and NLO QCD calculations obtained with 
the program MEPJET give a good description of the measurements of inclusive and dijet 
cross sections. 

The average number of subjets decreases as E 3 ^ 1 increases. The NLO QCD calculations 
agree well with the measured subjet multiplicities, (ri s bj)- A fit of the measured (n s bj) as 
a function of E 3 ^ at y cut = 10 -2 provides a determination of the strong coupling constant 
a s (Mz). The value of a s (Mz) determined for the region Ef > 25 GeV is 

a s (M z ) = 0.1202 ± 0.0052 (stat.)iS (syst.) ±Hol (th.)- 
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This result is consistent with other recent determinations and with the PDG value. 

The subjet multiplicities in CC and NC DIS are similar as a function of Q 2 . The measured 
(n s bj) at a given is somewhat smaller in CC DIS than in NC DIS. This can be 
attributed to the different Q 2 distributions of the two processes. 
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1995-1997 e+p data 


sample f-v/s 


= 300 GeV) 




Q 2 range 
(GeV 2 ) 


da/dQ 2 
(pb/GeV 2 ) 


Astat 


^syst 




QED 
correction 


Chad 


200 - 500 


0.0105 


±0.0010 


+0.0012 
-0.0010 


+0.0014 
-0.0011 


1.045 


0.985 


500 - 1000 


0.0108 


±0.0008 


+0.0007 
-0.0005 


+0.0004 
-0.0004 


1.033 


1.000 


1000 - 2000 


0.00571 


±0.00042 


+0.00014 
-0.00037 


+0.00006 
-0.00003 


1.042 


1.001 


2000 - 4000 


0.00221 


±0.00018 


+0.00005 
-0.00008 


+0.00007 
-0.00004 


1.064 


0.999 


4000 - 10000 


0.000380 


±4.1 • 10" 5 


+1.0 . 1Q -5 


+3.6 i n -5 
-3.0 ' iu 


1.085 


0.998 





1999-2000 e + p data 


sample 


= 318 GeV) 




Q 2 range 
(GeV 2 ) 


da/dQ 2 
(pb/GeV 2 ) 


A s tat 


^syst 




QED 

correction 


Chad 


200 - 500 


0.0125 


±0.0010 


+0.0012 
-0.0014 


+0.0016 
-0.0013 


1.055 


0.983 


500 - 1000 


0.0107 


±0.0007 


+0.0007 
-0.0006 


+0.0005 
-0.0005 


1.048 


0.999 


1000 - 2000 


0.00668 


±0.00038 


+0.00010 
-0.00014 


+0.00009 
-0.00007 


1.054 


1.000 


2000 - 4000 


0.00233 


±0.00016 


+0.00002 
-0.00008 


+0.00006 
-0.00004 


1.073 


0.998 


4000 - 10000 


0.000489 


±4.0 ■ 10~ 5 


+0.9 ln -5 
-1.3 ' iu 


+ 4i ■ io 5 


1.097 


0.997 



Combined 1995-2000 e + p data 


sample (a/s 


= 318 GeV) 


Q 2 range 


da/dQ 2 








CT 3i8/°300 


(GeV 2 ) 


(pb/GeV 2 ) 


A s tat 


^syst 




200 - 500 


0.0119 


±0.0007 


+0.0012 
-0.0012 


+0.0016 
-0.0013 


1.0476 


500 - 1000 


0.0110 


±0.0005 


+0.0007 
-0.0005 


+0.0005 
-0.0004 


1.0580 


1000 - 2000 


0.00647 


±0.00029 


+0.00008 
-0.00018 


+0.00008 
-0.00005 


1.0750 


2000 - 4000 


0.00237 


±0.00012 


+0.00002 
-0.00007 


+0.00007 
-0.00005 


1.1022 


4000 - 10000 


0.000472 


±3.1 • 10^ 5 


+0.7 . in-5 
-1.2 iu 


+4.2 . , n -5 
-3.4 iu 


1.1633 



Table 1: Inclusive jet cross-section da/dQ 2 for jets of hadrons in the laboratory 
frame. The statistical, systematic and energy-scale uncertainties are shown sepa- 
rately. The multiplicative correction applied to correct for QED radiative effects 
and for hadronisation effects and the theoretical correction factor used to combine 
the two data sets are shown. 
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1995-1997 e + p data sample (y/s = 300 GeV) 


rf et range 


da / dff et 
(pb) A stat A sys t A ss 


QED 

correction 


Chad 


-1 toO 

to 1 

1 to 1.5 
1.5 to 2 


3.80 ±0.37 ±gS t°oit 

9.60 ±0.53 

inM -|-n 7Q +0-22 +0.05 

ru.oo rnu./y _ 38 _ 004 
7.78 ±0.67 ±8^ 


1.082 
1.052 
1.045 
1.042 


0.967 
0.992 
1.005 
1.017 




1999-2000 e + p data sample (y/s = 318 GeV) 


rf et range 


(pb) A st at A syst A^5 


QED 

correction 


Chad 


-1 toO 

to 1 

1 to 1.5 
1.5 to 2 


4 7fi + n otr +0.11 +0.20 

<±. zru.oo _ 044 _ 015 

q 7^ -+-n +°- 38 + - 06 
y.io _ 23 _ 004 

n.46 ±0.69 ±8:81 

n ct; J_n fiO +0.25 +0.04 
y.jo zcu.uz _ 029 _ _ 04 


1.086 
1.064 
1.058 
1.060 


0.967 
0.991 
1.003 
1.013 



Combined 1995-2000 e + p data sample (y/s = 318 GeV) 


rf ct range 


da/d7f et 
(pb) A stat A syst A ES 


^s/^ioo 


-1 toO 

to 1 

1 to 1.5 
1.5 to 2 


4.48 ±0.26 ±8^g ±8:ii 

9.98 ±0.35 ig : « 

ii.so ±0.54 Tl ±8:8S 

9.26 ±0.48 ±«:|? ±g : gg 


1.0670 
1.0797 
1.0990 
1.1299 



Table 2: Inclusive jet cross-section da/dyf et for jets of hadrons in the laboratory 
frame. For details, see the caption to Tabled 
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1995-1997 e + p data sample 


(x/i = 300 GeV) 




range 


dajdEf 








QED 


r 

^had 


(GeV) 


(pb/GeV) 


A 

^stat 


<—^syst 




correction 


14 - 21 


0.950 


±0.068 


+0.081 
-0.069 


+0.014 
-0.014 


1.025 


0.992 


21 - 29 


0.703 


±0.052 


+0.040 
-0.051 


+0.006 
-0.004 


1.036 


0.998 


29 - 41 


0.486 


±0.034 


+0.008 
-0.022 


+0.006 
-0.004 


1.066 


0.998 


41 - 55 


0.219 


±0.021 


+0.008 
-0.007 


+0.004 
-0.004 


1.094 


0.998 


55 - 71 


0.0542 


±0.0096 


+0.0026 
-0.0026 


+0.0020 
-0.0020 


1.131 


0.988 


71 - 87 


0.0210 


±0.0058 


+0.0007 
-0.0008 


+0.0011 
-0.0012 


1.147 


0.986 


87 - 119 


0.00643 


±0.00214 


+0.00121 
-0.00064 


+0.00070 
-0.00047 


1.219 


0.981 





1999-2000 


e + p data sample 




GeV) 




E^ range 
(GeV) 


dajdEf 
(pb/GeV) 


^stat 


^syst 




QED 

correction 


Chad 


14 - 21 


1.030 


±0.060 


+0.056 
-0.092 


+0.014 
-0.012 


1.035 


0.989 


21 - 29 


0.816 


±0.047 


+0.022 
-0.034 


+0.008 
-0.007 


1.052 


0.995 


29 - 41 


0.527 


±0.030 


+0.014 
-0.007 


+0.006 
-0.005 


1.071 


0.998 


41 - 55 


0.230 


±0.018 


+0.003 
-0.006 


+0.004 
-0.004 


1.112 


0.998 


55 - 71 


0.0775 


±0.0096 


+0.0026 
-0.0035 


+0.0026 
-0.0026 


1.135 


0.990 


71 - 87 


0.0232 


±0.0052 


+0.0010 
-0.0012 


+0.0012 
-0.0012 


1.172 


0.986 


87 - 119 


0.00385 


±0.00146 


+0.00034 
-0.00012 


+0.00036 
-0.00033 


1.192 


0.984 



Combined 1995-2000 e + p data 


sample 


(Vs = 318 


GeV) 


£y range 


dajdEf 








a 318/°300 


(GeV) 


(pb/GeV) 


Astat 


^syst 




14 - 21 


1.022 


±0.046 


+0.067 
-0.084 


+0.014 
-0.013 


1.0621 


21 - 29 


0.791 


±0.036 


+0.027 
-0.034 


+0.007 
-0.006 


1.0688 


29 - 41 


0.527 


±0.023 


+0.009 
-0.009 


+0.006 
-0.005 


1.0855 


41 - 55 


0.236 


±0.014 


+0.004 
-0.004 


+0.004 
-0.004 


1.1216 


55 - 71 


0.0727 


±0.0074 


+0.0020 
-0.0029 


+0.0025 
-0.0025 


1.1891 


71 - 87 


0.0246 


±0.0043 


+0.0010 
-0.0008 


+0.0013 
-0.0014 


1.3046 


87 - 119 


0.00572 


±0.00143 


+0.00066 
-0.00027 


+0.00058 
-0.00045 


1.5461 



Table 3: Inclusive jet cross-section da/dE^ for jets of hadrons in the laboratory 
frame. For details, see the caption to Tabled 
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1995-1997 e + p data 


sample (y^s 


= 300 GeV) 






Q 2 range 
(GeV 2 ) 


da/dQ 2 
(pb/GeV 2 ) 




^syst 




QED 

correction 


Chad 


200 - 500 
500 - 1000 
1000 - 2000 


0.00153 
0.00165 
0.00112 


±0.00051 
±0.00034 
±0.00020 


+0.00030 
-0.00067 
+0.00015 
-0.00022 
+0.00006 
-0.00007 


+0.00031 
-0.00023 
+0.00014 
-0.00013 
+0.00003 
-0.00003 


1.040 
1.045 
1.059 


0.916 
0.924 
0.926 


2000 - 4000 


3.85 • 10" 4 


±0.82 • 10" 4 


+0.25 I n -4 
-0.26 ' iu 


+0.26 1 n -4 
-0.15 ' iu 


1.077 


0.910 


4000 - 10000 


7.13 • 10" 5 


±1.91 • 10" 5 


+0.70 I n -5 
-0.12 ' iu 


+0.77 n -5 
-0.79 ' iu 


1.083 


0.893 





1999-2000 e + p data sample (y/s = 


318 GeV) 






Q 2 range 
(GeV 2 ) 


da/dQ 2 
(pb/GeV 2 ) 


A s tat 


^syst 




QED 

correction 


Chad 


200 - 500 


0.00290 


±0.00062 


+0.00081 
-0.00079 


+0.00057 
-0.00043 


1.068 


0.914 


500 - 1000 


0.00190 


±0.00031 


+0.00023 
-0.00026 


+0.00017 
-0.00015 


1.054 


0.925 


1000 - 2000 


0.00112 


±0.00016 


+0.00006 
-0.00007 


+0.00004 
-0.00003 


1.067 


0.927 


2000 - 4000 


4.02 • 10~ 4 


±0.70 ■ 10~ 4 


+0.27 1fl -4 
-0.32 ' iu 


+0.16 i n -4 
-0.15 ' iu 


1.078 


0.916 


4000 - 10000 


11.70 • 10~ 5 


±2.04 • 10" 5 


+1.12 1fl -5 
-0.99 ' iu 


+1.32 ln -5 
-0.99 ' iu 


1.101 


0.904 



Combined 1995-2000 e + p data sample (y/s = 318 GeV) 


Q 2 range 
(GeV 2 ) 


da/dQ 2 

(pb/GeV 2 ) A siai A syst A ES 


th / th 

°318/ °300 


200 - 500 
500 - 1000 
1000 - 2000 
2000 - 4000 
4000 - 10000 


0.00242 ±0.00043 1™ 

0.00186 ±0.00024 t° Zll 

00116 ±0 00013 +0.00006 +0.00004 

U.UUIIU ICU.UUUIO -0.00005 -0.00003 

4.i3 -lo- 4 ±0.56 -lo- 4 ±8;i| - 10- 4 ±o: 2 6-io- 4 

10.49 • 10" 5 ±1.53 ■ 10" 5 t° f 7 ■ 10" 5 toll ■ 10 " 5 


1.0848 
1.0859 
1.0926 
1.1191 
1.1777 



Table 4: Dijet cross-section da/dQ 2 for jets of hadrons in the laboratory frame. 
For details, see the caption to Tabled 
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1995-1997 e + p data sample (y/s = 300 GeV) 


mi2 range 
(GeV) 


da/dm-12 

(pb/GeV) A stat A syst A ES 


QED 
correction 


Chad 


10 - 20 
20 - 30 
30 - 40 
40 - 75 


0.0724 ±0.0137 j£882S + -°oZl 

o.io4 ±0.018 i^ 6 7 ±8:881 

0.0884 ±0.0184 ^ %™ 
0.0296 ±0.0068 ±8l8gS 


1.050 
1.053 
1.060 
1.089 


0.921 

0.901 
0.918 
0.942 




1999-2000 e + p data sample (y/s = 318 GeV) 


m 12 range 
(GeV) 


da /dmi2 

(pb/GeV) ^stat A syst A ES 


QED 
correction 


Chad 


10 - 20 
20 - 30 
30 - 40 
40 - 75 


0.0913 ±0.0129 t° S + -°oS 

o.i4o ±0.018 ™» ±8:8St 
o.no ±0.018 t° ™ ±8:881 

0.0228 ±0.0049 ±8'oo2i ^oon 


1.059 
1.064 
1.082 
1.085 


0.918 
0.902 
0.917 
0.941 



Combined 1995-2000 


e + p data 


sample 


(y/s = 318 GeV) 


77ii2 range 


da jdmyi 








th 1 th 
°318/ °300 


(GeV) 


(pb/GeV) 


^■stat 


^syst 




10 - 20 


0.0863 


±0.0098 


+0.0058 
-0.0057 


+0.0018 
-0.0016 


1.0830 


20 - 30 


0.130 


±0.013 


+0.006 
-0.009 


+0.003 
-0.003 


1.0992 


30 - 40 


0.106 


±0.013 


+0.009 
-0.008 


+0.005 
-0.003 


1.1172 


40 - 75 


0.0270 


±0.0042 


+0.0051 
-0.0031 


+0.0013 
-0.0013 


1.1552 



Table 5: Dijet cross-section da/dm^ for jets of hadrons in the laboratory frame. 
For details, see the caption to Tabled 
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y cut value 


(n sbj ) 


A 


stat 


^syst 


Chad 


14 < < 17 GeV 


0.0005 


4.432 


±0 


082 


+0.078 
-0.025 


2.127 


0.001 


3.576 


±0 


071 


+0.098 
-0.020 


1.845 


0.003 


2.616 


±0 


057 


+0.058 
-0.020 


1.526 


0.005 
0.01 


2.238 
1.755 


±0 
±0 


050 
042 


+0.039 
-0.022 

+0.026 
-0.031 


1.426 
1.292 


0.03 


1.263 


±0 


033 


+0.028 
-0.015 


1.094 


0.05 


1.117 


±0 


023 


+0.019 
-0.013 


1.037 


n 1 

U. 1 


i nm 


±0 


007 


+0.018 
-0.003 


l.UUO 




17 < 


< 21 GeV 




0.0005 


4.272 


±0 


071 


+0.067 
-0.050 


1.989 


0.001 


3.522 


±0 


060 


+0.058 
-0.043 


1.726 


0.003 
0.005 


2.427 
1.999 


±0 
±0 


047 
042 


+0.049 
-0.025 

+0.036 
-0.012 


1.464 
1.364 


0.01 


1.574 


±0 


034 


+0.021 
-0.017 


1.227 


0.03 


1.176 


±0 


026 


+0.014 
-0.018 


1.054 


0.05 
0.1 


1.102 
1.021 


±0 
±0 


020 
009 


+0.021 
-0.010 

+0.006 
-0.006 


1.016 
0.999 



Table 6: Mean subjet multiplicity as a function of y cnt for the E^ regions 14 < 
El£ < 17 and 17 < E^ < 21 GeV. The statistical and systematic uncertainties are 
shown separately. The multiplicative correction applied to correct for hadronisation 
effects is shown in the last column. 
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y cut value 


(n sbj ) 


A 


stat 


^syst 


Chad 




21 < 


< 25 GeV 




U.UUUo 


4 


051 


±0 


070 


+0.067 
-0.010 


l.OO 1 


0.001 


3 


251 


±0 


059 


+0.071 
-0.011 


1.637 


0.003 


2 


266 


±0 


043 


+0.052 
-0.014 


1.409 


0.005 


1 


953 


±0 


038 


+0.042 
-0.008 


1.304 


0.01 


1 


504 


±0 


034 


+0.035 
-0.007 


1.167 


0.03 


1 


164 


±0 


025 


+0.023 
-0.008 


1.025 


0.05 


1 


070 


±0 


017 


+0.023 
-0.005 


1.003 


0.1 


1 


022 


±0 


009 


+0.008 
-0.005 


0.997 


25 < < 35 GeV 


0.0005 


3 


786 


±0 


049 


+0.027 
-0.021 


1.736 


0.001 


2 


997 


±0 


039 


+0.028 
-0.011 


1.570 


0.003 


2 


067 


±0 


029 


+0.017 
-0.009 


1.339 


0.005 


1 


717 


±0 


025 


+0.019 
-0.012 


1.228 


0.01 


1 


386 


±0 


023 


+0.018 
-0.009 


1.099 


0.03 


1 


137 


±0 


016 


+0.005 
-0.012 


1.006 


0.05 


1 


061 


±0 


Oil 


+0.003 
-0.013 


0.998 


0.1 


1 


008 


±0 


004 


+0.004 
-0.004 


0.998 



Table 7: Mean subjet multiplicity as a function of y cut for the E^ regions 
21 < Ej? < 25 and 25 < E^ < 35 GeV. For details, see the caption to Table® 
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y cut value 


(n sbj ) 


A 


stat 


^syst 


Chad 




35 < 4 ct 


< 55 GeV 




U.UUUO 


3 


343 


±0 


045 


+0.037 
-0.017 


l.o < y 


0.001 


2 


649 


±0 


040 


+0.038 
-0.013 


1.433 


0.003 


1 


797 


±0 


029 


+0.029 
-0.009 


1.197 


0.005 


1 


512 


±0 


026 


+0.027 
-0.006 


1.106 


0.01 


1 


245 


±0 


022 


+0.020 
-0.004 


1.026 


0.03 


1 


106 


±0 


015 


+0.003 
-0.002 


0.996 


0.05 


1 


060 


±0 


Oil 


+0.005 
-0.001 


0.996 


0.1 


1 


010 


±0 


004 


+0.001 
-0.000 


0.998 


55 < E^ < 119 GeV 


0.0005 


2 


790 


±0 


068 


+0.013 
-0.035 


1.450 


0.001 


2 


196 


±0 


056 


+0.022 
-0.026 


1.285 


0.003 


1 


520 


±0 


041 


+0.038 
-0.012 


1.073 


0.005 


1 


378 


±0 


042 


+0.040 
-0.015 


1.022 


0.01 


1 


260 


±0 


039 


+0.024 
-0.012 


0.998 


0.03 


1 


098 


±0 


025 


+0.005 
-0.007 


0.996 


0.05 


1 


044 


±0 


017 


+0.005 
-0.003 


0.998 


0.1 


1 


036 


±0 


026 


+0.004 
-0.002 


1.000 



Table 8: Mean sub jet multiplicity as a function of y cut for the regions 
35 < Ej? < 55 and 55 < E^ < 119 GeV. For details, see the caption to Table® 
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Q 2 range 
(GeV 2 ) 


(n sbj ) 


Astat 


^syst 


Chad 


200 - 350 


1.695 


±0.053 


+0.025 
-0.027 


1.316 


350 - 500 


1.677 


±0.050 


+0.011 
-0.049 


1.257 


500 - 750 


1.516 


±0.036 


+0.028 
-0.011 


1.212 


750 - 1000 


1.487 


±0.040 


+0.042 
-0.018 


1.175 


1000 - 2000 


1.391 


±0.024 


+0.030 
-0.012 


1.129 


2000 - 4000 


1.372 


±0.027 


+0.012 
-0.014 


1.081 


4000 - 10000 


1.318 


±0.031 


+0.034 
-0.010 


1.051 



Table 9: Measurement of the mean subjet multiplicity at y cnt = 10 2 as a function 
of Q 2 . For details, see the caption to Table® 
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Figure 1: (a) The differential cross-section da/dQ 2 for inclusive jet production 
in the laboratory frame with > 14 GeV and —1 < rf et < 2 in the kinematic 
region Q 2 > 200 GeV 2 and y < 0.9 for the 1995-2000 e + p data (black dots). The 
data are corrected to hadron level. The inner error bars represent the statistical 
uncertainty of the data, the outer error bars show the statistical and the system- 
atic uncertainties (not associated with the uncertainty in the absolute energy scale) 
added in quadrature. The shaded band displays the uncertainty due to the abso- 
lute energy scale of the GAL. The parton shower Monte Carlo prediction given by 
ARIADNE at hadron level (dashed line) and the next-to-leading- order prediction 
obtained with MEPJET corrected to hadron level (solid line) are shown, (b) The 
ratio of the measured cross section to the next-to-leading- order calculation. The 
theoretical uncertainty is indicated by the hatched band. 
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Figure 2: (a) The differential cross-section da /drf et for inclusive jet production 
in the laboratory frame with > 14 GeV and — 1 < rf ei < 2 in the kinematic 
region Q 2 > 200 GeV 2 and y < 0.9. Other details are as decribed in the caption 
to Fig. □ 
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Figure 3: (a) The differential cross-section da/dEjT 1 for inclusive jet production 
in the laboratory frame with E^ 1 > 14 GeV and — 1 < rf ei < 2 in the kinematic 
region Q 2 > 200 GeV 2 and y < 0.9. Other details are as decribed in the caption 
to Fig. □ 
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Figure 4: (a) The differential cross-section da/dQ 2 for dijet production in the 
laboratory frame with E^' > 14 GeV, E 3 ^ 1 ' 2 > 5 GeV and —1 < rf et < 2 in the 
kinematic region Q 2 > 200 GeV 2 and y < 0.9. Other details are as decribed in the 
caption to Fig. 
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Figure 5: (a) The differential cross-section da/dmi2 for dijet production in the 
laboratory frame with E^ 1 ' 1 > 14 GeV, E 3 ^ 1 ' 2 > 5 GeV and — 1 < rj jet < 2 in the 
kinematic region Q 2 > 200 GeV 2 and y < 0.9. Other details are as decribed in the 
caption to Fig. [7J 
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Figure 6: Mean subjet multiplicity (black dots), (n s bj) , as a function of y cut for 
inclusive jet production in the laboratory frame with — 1 < rj^ 1 < 2 in different E^ 1 
regions. The parton shower Monte Carlo prediction given by ARIADNE at hadron 
level (dashed line) and the next-to-leading- order prediction obtained with MEPJET 
corrected to hadron level (solid line) are shown. The error bars are smaller than 
the symbols. 
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Figure 7: Mean subjet multiplicity, (n s bj), at y cut = 10 2 as a function of E 3 ^ 
(black dots), for inclusive jet production in the laboratory frame with > 14 
GeV and — 1 < rf et < 2. The inner error bars represent the statistical uncertainty 
of the data. The outer error bars show the statistical and systematic uncertainties 
added in quadrature. The NLO QCD predictions obtained with MEPJET using the 
CTEQ4 sets of proton PDFs are shown for 3 different values of a s (M z ) (curves). 
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Figure 8: Measurements of (n s bj) at y cu t = 10 -2 for inclusive jet production in 
the laboratory frame with > 14 GeV and — 1 < rj jet < 2 in charged current 
DIS (circles) and neutral current DIS (open squares) as a function of (a) E 3 ^ and 
(b) Q\ 



32 



